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^ ■ Abstract 



In a recent publication we proposed an extension of Hirota's bilinear formalism to arbitrary 
multilinearities. The trilinear (and higher) operators were constructed from the requirement 
of gauge invariance for the nonlinear equation. Here we concentrate on the trilinear case, 
and use singularity analysis in order to single out equations that are likely to be integrable. 
New PDE's are thus obtained, along with others already well-known for their integrability 
and for which we obtain here the trilinear expression. 
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1. Introduction 

What is astonishing with integrable PDE's is that there are so many of them. Infi- 
nite hierarchies of equations are known to date and the domain is still expanding. 
Moreover, these equations are often associated to physical models [1]. Then, us- 
ing integrability, one can construct particular solutions and compute the pertinent 
physical quantities for these models. 

One of the essential tools in the study of integrable PDE's, over the past 20 
years, has been the bilinear formalism of Hirota [2]. The main advantage of this 
approach lies in the fact that it allows one to obtain multisoliton solutions in a 
straightforward way [3] . Still, the proof of integrability, based on the existence of 
an arbitrary number of solitons, may present considerable difficulties. It is often 
easier to test the integrability of a given equation using singularity analysis (Painleve 
method) [4] , especially since this criterion (absence of singular solutions that exhibit 
branching) can be implemented algorithmically. 

In a recent work [5] we have presented an extension of Hirota's bilinear formalism 
to higher multilinearities. This was motivated by the existence of integrable equa- 
tions that cannot be bilinearized, like the Satsuma determinantal trilinear PDE's 
[6] . (It must be made clear at this point that the nonexistence of a bilinear form is 
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not a rigorous statement; it simply means that no bilinear form has been obtained 
to date). Our method is based on the assumption of gauge invariance. In the first 
exploratory study [5] we studied only the leading part of one-dimensional trilinear 
equations with one dependent variable and found those that satisfy the Painleve 
criterion. The aim of the present paper is to complete this work, and obtain the 
integrable higher dimensional generalizations with the non-leading parts. As we 
will show in the following sections some of the equations obtained are new, while 
others are well-known for their integrability without necessarily possessing a simple 
bilinear form. 

2. From bilinear to trilinear operators. 

A prerequisite to the application of the Hirota bilinear formalism is a dependent 
variable transformation that converts the nonlinear equation into a quadratic 'pre- 
potential' form. In order to make things more clear let us present the classical 
example of the KdV equation. Starting from 

u xxx + Quu x + u t = 0, (2.1) 

we introduce the transformation u = 2d x logF and obtain (after one integration): 

F XXXX F — AF XXX F X + 3F XX + F xt F — F x F t = 0. (2.2) 

This last equation can be written in a particularly condensed form using the Hirota 
D operator: 

(D x + D x D t )F-F = 0, (2.3) 

where the D-operator is defined by its antisymmetric derivative action on a pair of 
functions (the 'dot product'): 

D n x f-9 = (d Xl -d xa ) n f(x 1 )g(x 2 )\ xa=xi=x . (2.4) 

The crucial observation here is the relation of the 'physical' variable u = 2d x log F 
to the Hirota's function F: the gauge transformation F — » e px+lAjt F leaves u invari- 
ant. It turns out that this is a general property of bilinear equations. In fact, one 
can define Hirota's bilinear equations through the requirement of gauge invariance. 
This statement was proven in [5]. 

Having obtained the Hirota bilinear operators on the sole requirement of gauge 
invariance we investigated in [5] the possible extension of this formalism and the 
introduction of multilinear operators. This turned out to be possible. Our first 
step was the extension of the bilinear operators to the trilinear case. We found 
that one convenient basis for the representation of the gauge invariant AT'th order 
derivative operators is given by (d\ — <92) n (<9i — ds) N ~ n for n = 0, . . . , N, where 
the subscripts tell on which of the three functions each derivative acts. Thus the 
basic building blocks for the trilinear operators are again the Hirota bilinear .D's: 



we must just specify the indices in this case. We thus have Di 2 = d Xl — d X2 , D 23 = 
d X2 ~~ 9 X3 , -D31 = d X3 — d Xl , but, of course the three are not linearly independent: 
D12 + -D23 + -D31 = 0. Their action on a 'triple dot product' is analogous to the 
bilinear case: 

D 12 f-g-h = (d Xl - d X2 )f(x 1 )g(x 2 )h(x 3 )\ x:i=x2=xi=x = (f'g - fg')h. (2.5) 

The choice of a particular pair of .D's as the basic trilinear operators breaks the 
symmetry between the three coordinates x^s. It is possible to restore this symmetry 
by introducing a different basis for the trilinear operators, T and T*: 

T = d 1 + jd 2 + fd 3 , T* = d x + fd 2 + jd 3 , (2.6) 

where j is the cubic root of unity, j = e 2 * 71 "/ 3 . (Note that the star in T* indicates 
complex conjugation for the coefficients in T but not for the independent variables). 
Note that T n T* m F-F-F = 0, unless n — m = (mod 3), which is the equivalent 
to the bilinear property D n F-F = 0, unless n = (mod 2). Moreover, a reality 
condition is satisfied: T n T* m F-F-F = T m T* n F-F-F. 

The generalization to higher multilinear equations is straightforward. One can 
introduce the set of n(n — l)/2 operators Dij acting on n-tuple dot-products Dijfi- 
fi- ■ ■ ■ -fn- Of course only n — 1 of the D^ 's are independent, a convenient basis being 
the -Dij, j — 2, . . .n. As in the trilinear case, one can also construct 'symmetric' 
operators: 

n-l 

M™ = J2 e 27Ttkm/n d k+1 (2.7) 

k=0 

for < m < n. We have, for example, D = M 2 , T = M3, T* — M% and so on. 

3. General comments on mult ilinearizat ion 

Before dealing with specific cases let us present here some general considerations. 
Let us start with a nonlinear (in u) equation, of order k + 1 having the form 
ut + d x P(u, u x , . . . , Ukx) = 0. Several well known integrable equations belong to 
this class. Let us assume that the leading part of P is weight-homogeneous in u 
and d x , with u having the same weight as d 2 . Then we can transform the equation 
into a multilinear expression through the transformation u = ad 2 (log F) and ob- 
tain, after one rr-integration, generically an (/c + 2)-multilinear equation. The scaling 
factor a can then be chosen (perhaps in several ways) so as to make an F 2 term 
factor out: the resulting multilinear equation is at most /c-linear. For example at 
order five (k = 4) we have three integrable equations, and we should expect, in prin- 
ciple, these equations to have quadrilinear forms. Some unexpected cancellations, 
however, do occur. Thus the Sawada-Kotera equation [7] has a bilinear expression, 
the Lax-5 equation [8] a trilinear form, but the Kaup-Kuperschmidt equation [9] is 
quadrilinear. 
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In the process of multilinearization of a given nonlinear equation the following 
points should be noted: 

1) Above we discussed only u = ad 2 log(F) substitutions. If u always appears in 
the equation with at least one derivative then we can get a gauge invariant form also 
with u = adlog(F), and if we always have at least two derivatives, u = alog(F) is 
sufficient. 

2) Any bilinear expression multiplied by F is cubic and gauge invariant and 
therefore has trilinear form, for example 

(3.1) 

(3.2) 

(3.3) 

)F-F-F. (3.4) 

Therefore known one-component bilinear equations reappear in this trilinear study. 

3) Taking a derivative of an equation having an n-linear form yields an equation 
that is n + 1 multilinear, for example 

d x (F~ 2 D x D y F-F) = |F- 3 T x 2 T y F-F-F, (3.5) 

d x (F- 2 D 3 x D y F-F) = |F" 3 (T*T; + 2T 3 T*T y ) F-F-F, (3.6) 

d x (F- 2 D 6 x F-F) =lF~ 3 TlTfF-F-F. (3.7) 

The integr ability of the derivative form is another matter: if for example the x- 
derivate form is also integrable then the original equation would be integrable with 
inhomogeneous ^-independent term. Thus we will later find the x-derivative of 
Ito equation to be integrable, but not its t-derivative, therefore Ito's equation is 
integrable even with a inhomogeneous g(t)-term. 

Since equations having a bilinear form will reappear in the trilinear analysis as 
mentioned before, let us collect here the integrable equations that can be bilinearized 
to the form P(D)F-F = [10,11]: 

Kadomtsev-Petviashvili: 

(D 4 x -4D x D t +3D 2 y )F-F = 0, (3.8) 

Ito: 

(D t D s x + aD 2 x + D t D y )F-F = 0, (3.9) 

Hietarinta: 

{D x (Dl - D 3 t ) + aD 2 x + bD t D x + cD 2 ) F-F = 0, (3.10) 

Sawada-Kotera-Ramani: 

(D 6 X + 5D 3 x D y - 5D 2 y + D x D t ) F-F = 0, (3.11) 



4. Singularity analysis of trilinear equations 

One motivation behind the multilinear approach is to find new integrable equations. 
Eventually we would like to get something similar to the systematic classification 
of bilinear equations presented in [10,11]. In the following paragraphs we will limit 
ourselves to the singularity analysis of trilinear equations involving only one depen- 
dent variable, i.e. unicomponent equations P(T x ,T*,Ty,T*, . . .)F-F-F = 0. (Let 
us recall here that in the bilinear case a complete classification of these simplest 
equations [10,11] was possible). 

Since the dependent function u in a nonlinear equation is related to the mul- 
tilinear F through u = ad 2 log F it is clear that a zero in F induces a pole-like 
behavior in u. Let us here consider a concrete example: T 3k F-F-F = 0. Putting, 
for the dominant part, F ~ <p n , where <p is the singular manifold we find that the 
only possible leading behaviors have n — 0, 1, . . .,k — 1. The first corresponds to a 
nonsingular Taylor-like expansion, which is always possible. The second behavior 
F ~ (/) corresponds to a simple zero of F that would give a (double) pole in u. 
Next the resonances r can be obtained if we substitute F = <p n (l +u)(j) r ) and collect 
terms linear in to. For the equation to pass the Painleve test the resonances must 
be integers and no incompatibilities must arise at any resonance. As an example 
let us take k = 2, i.e., the operator T 6 . The leading behavior F ~ leads to the 
resonances r = —1, 0, 1, 2 and the roots of r 2 — 13r + 60 = which are complex. 
Thus T 6 F-F-F = does not possess the Painleve property. 

4.1 Leading behaviors 

The Painleve analysis for the leading part of the equations was performed in [5] , with 
just one independent variable. This is sufficient for the computation of dominant 
singularities and resonances. Note that for a given N = n + m there may exist 
several pairs (n, m) such that T n T* m F-F-F is not identically zero, namely those 
for which n = m (mod 3). The leading part of the general equation at order N 
is then given by a linear combination of all the non- vanishing T n T* rn F-F-F , s. In 
each case we give below the precise combinations that lead to equations with the 
Painleve property. The nonlinear forms of the leading parts of the equations are 
obtained by the standard substitution F = e 9 followed by u = 2g". The results of 
[5] are summarized as follows: 
N = 2 : operator : TT* leading part : u. 

In this case we can write the result also in bilinear form 
N = 3 : operator : T 3 leading part : u'. 
N = 4 : operator : T 2 T* 2 leading part : u" + 3u 2 . 

This, of course, is just the leading part of the KdV equation in potential form, 

and can also be written in bilinear form, see (3.2). 
N = 5 : operator : T 4 T* leading part : v!" + 6u'u. 

This is the derivative of the expression obtained at iV = 4. 
N = 6 : operator : AT 6 + / uT 3 T* 3 . This is the first case where we have two possible 
(n, m) pairs. The A, [i combinations that pass the Painleve test are the following 



a) (A = 7, fj, = 20) leading part : u"" + 10u"u + bu' 2 + 10u 3 . 

This is the leading part of the once integrated 5th order equation in the Lax 
hierarchy of KdV. 

b) (A = -2,/jl= 20) leading part : u"" + 15u"u + 15u 3 . 

This is the leading part of the once integrated Sawada-Kotera-Ramani equation 
(3.11), and can also be written in bilinear form, see (3.3). 

c) (A = —/j,) leading part : uu" — u' 2 + u 3 . 

This expression can be cast in determinantal equation 

pim pin pi/ 
pin pn pi 

F" F' F 

whose integrability was already noticed by Chazy [12]. 
N = 7: operator : T 5 T* 2 leading part : w (5) + 15u"'u + 15u"u' + 45u'u 2 . 

This is the Sawada-Kotera equation, i.e., the derivative of the expression obtained 

iniV = 6b, c.f. (3.7). 
N = 8 : operator : XT 7 T* + (iT^T* 4 and nonlinearization with u = 6g" instead 
of u = 2g" used before. 

a) (A = 4, fj, = 5) leading part : u^+6u""u+10u'"u'+5u" 2 + 10u"u 2 +10u' 2 u+^u 4 
This would correspond to a 7 order equation which, we believe, leads to a new 
integrable case. 

b) (A = 4,ji= 14) leading part : u^ +7u""u+7u'"u' + 7u" 2 + Uu"u 2 + 7u' 2 u+lu 4 
This is the leading part of a once integrated higher Sawada-Kotera equation. 

c) (A = -/z) leading part : u""u - 3u'"u' + 2u" 2 + 4u"u 2 - 3u ,2 u + \u A . 

This new equation looks like an extension of Satsuma's equation given at iV=6c 

above, but it cannot be written as a single determinant. 
N = 9: No cases passing the Painleve test exist at this order. 
N = 10 : operator : 5T 8 T* 2 + 4T 5 T* 5 (and we take u = 30g"). 

leading part : u^ +2u^u+Au^u' + Qu""u" + hu"' 2 + \ u""u 2 +Au'"u'u+2u" 2 u + 

2u"u' 2 + ±u"u 3 + \u" 2 u 2 + ^w 5 

This is also a new equation. 
Furthermore, no integrable candidates were found at orders iV = 11, 12. The singu- 
larity analysis at these higher orders becomes progressively more difficulty Already 
at iV = 12 there exist three non- vanishing n, m combinations). It is, thus, not 
possible to extend our investigation to very high orders (as was done in our study 
of bilinear equations [11]), but we do believe that no further integrable candidates 
exist at higher orders. 

4.2 Resonance conditions 

In order to investigate the resonance conditions we must specify precisely the PDE 
we are working with, i.e., also the nonleading parts. Let us state from the outset 
that we are not going to consider equations of the form N = 6c or N = 8c. These 
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equations have a form that makes the expected leading singularity to vanish identi- 
cally, which leads to certain complications in the singularity analysis: they will be 
discussed elsewhere. Also in every case examined, subleading parts of opposite par- 
ity (= odd vs. even number of T"s) to that of the leading were found to violate the 
Painleve property, but we did not make an exhaustive study of this. The technical 
details of the resonance conditions study are not particularly interesting, the anal- 
ysis was performed on a computer using the REDUCE language [13] for symbolic 
manipulation. Although in some cases the resonances are particularly high (in the 
case of N = 10 we have a resonance at r = 30 for the leading behavior F ~ <p 14: ) we 
were able in all cases to check the resonance conditions. 

For low N the complete results are simple: for N = 3 one obtains only linear 
equations, and for N = 4 all the integrable equations could also be written in 
bilinear form. 

For equations with N > 5 an important question was how to extend the leading- 
order trilinear operator to two dimensions and still keep the calculations manage- 
able. Our choice has been to consider only 'monomial' leading parts, that is mono- 
mial in derivatives and not necessarily in T's. (From our experience on the bilinear 
case we do not expect nonmonomial leading parts to play a role for N > 5). 
N = 5 : 

The leading 1-dimensional operator is T^T* and three monomial generalizations 
to 1+1 dimensions are possible: i) T^T*, ii) T%T* + XT^T y T* and iii) T^T y T* + 

\rp2rp2rp* 

A x y x • 

The case i) supplemented by subleading terms leads only to the ^-derivative of 
the KP equation (3.5-8): 

(T*T* - 4T x 2 T t + 3T X T 2 )F-F-F = 0. 

In the case ii) (discarding the A = — 1 case that does not have a standard leading 
singularity) there exist the three Painleve candidates A = —4, 2 and 8. The first case 
does not satisfy the resonance condition. The case A = 2 yields just the ^-derivative 
of the Ito equation (3.9): 

(T^T* + 2T*T*T t + 3aT| + 3T 2 T y )F-F-F = 0. 

The final case A = 8 is more interesting. The most general result passing the 
Painleve test reads: 

(T y (T* 4 + 8T X 3 T X * + 9T y 2 ) + 9T x 2 T t ) F-F-F = 0. (4.1) 

Putting F = e 9 ' 2 we obtain: 

9xxxxy I ^9xxy9xx ' ^yxydxxx T 9yyy ' 9xxt "• V*'^) 

This equation is a generalization (the g yyy term is new) of an equation obtained by 
Schiff [14] from a reduction of the self-dual Yang-Mills equations. Incidentally, for 



this equation we have checked that a nontrivial three soliton solution does exist and 
this is a further indication of its integrability. It is interesting to observe that this 
equation reduces to the derivative KP of Case i) if the first factor T y is replaced by 
T x in (4.1). 

Finally the case iii) did not lead to any equations with the Painleve property. 

A general study of monomial leading parts was performed for all iV's higher than 
5. To make a long story short, only a monomial leading part of the form T^~ 1 T y 
has acceptable resonances for N = 6, 7. However, even in these cases, the resonance 
conditions are not satisfied. Thus, beyond N = 5 the leading part is not only 
monomial but also 1-dimensional. 
N = 6: 

The leading operator is XT® + 20T 3 T* with A =7 or —2. For A = 7 the equation 
obtained is the fifth-order equation in the Lax hierarchy [8] plus the KdV: 

(20T 3 T* 3 + ITl + aT 2 T* 2 + 27T*T y )F-F-F = (4.3) 

or in nonlinear form after one integration (u = 2d 2 logF): 

u 5x + 10uu xxx + 20u x u xx + 30u 2 u x + a(u xxx + 6uu x ) + u y = (4.4) 

Thus the Lax-5 equation does not possess a simple bilinear form like KdV but a 
trilinear one. For A = —2 the resulting equation is the 3-dimensional Sawada- 
Kotera-Ramani equation, 

(lOT^T* 3 -T%+ 45a(T 2 T*T* - aT y T*) + 9T X T Z )F-F-F = 0, (4.5) 

which also has a bilinear form (3.11). 

In the case N = 7 the only equation satisfying the Painleve requirement is the 
x-derivative of the Sawada-Kotera-Ramani equation, 

(3T^T* 2 + 5a(T^T y * - T 3 T*T y ) - 15a 2 T x T 2 - 3T 2 T Z )F-F-F = 0, (4.6) 

which means that the original equation is probably integrable even with an inho- 
mogeneous k(y, t) term. 

For N = 8 only the case 8a admits additive terms. Moreover it turns out that 
the first subleading term must also be 1-dimensional and it coincides with 6a. The 
full equation in this case reads: 

UT 7 X T* + 5T*T* 4 + a(20T 3 Tf + 7T X 6 ) + 9/3T 2 T*T* + ^-T X T*\ F-F-F = 0, 

(4.7) 
and its nonlinear form with u = 6d x log F is: 

u 7x + 6u 5x u x + 10u 4x u xx + hu 2 xxx + 10u xxx ul + -u x 

+a{3u 5x + 10u xxx u x + hu xx + —u x ) 

a/3 
+(3(Uxxy + u x u y ) + —u y = (4.8) 



To our knowledge this equation is a new integrable PDE. 

Case N = 10: Again one equation satisfies the integrability criterion. It reads: 

(5T|T; 2 + 4T X 5 T X * 5 + a(4T x 7 T x * + 5T X 4 T X * 4 ) + /3(20T X 3 T X * 2 T; + 7T x 5 T y ) 

+6a/3T^T*T* + -^-T y T*)F-F-F = (4.9) 

This equation, too, has not been encountered before as an integrable equation. 

5. Conclusion 

In the preceding paragraphs we have presented an extension of Hirota's bilinear 
formalism that can encompass any degree of multilinearity. The main guide in our 
investigation has been the requirement that the equations be gauge-invariant. Since 
our objective is the study of integrability we have also presented a classification of 
one-component trilinear equations that pass the Painleve test. 

An interesting difference between the bilinear and the tri- (and multilinear 
cases is that now free parameters enter already at the leading part. In practice 
this means that the Painleve analysis of the higher order unicomponent equations 
becomes increasingly difficult. Once the leading parts of these equations are fixed, 
we have studied the lower-order terms that can be added without destroying the 
Painleve property. As a result new integrable equations were discovered. The fact 
that they have eluded discovery till now is understandable since these equations are 
of high order and there is no hope to discover them by chance: a solid guide and a 
systematic approach are needed. 

Starting from a complete classification of unicomponent equations one can build 
up multicomponent ones following the approach we presented in [10,11] for the 
bilinear case. Further extensions can be presented and we can, of course, construct 
also higher multilinear equations (quadri-, penta-, etc.) equations. 

Extension of our formalism to the discrete case is also possible. In this case the 
integrability requirement is just the property of singularity confinement [15]. Some 
preliminary results exist already in this direction. 

Thus the multilinear extension to Hirota bilinear approach has a wide range of 
applicability. There are still many open problems, e.g., the computation of the mul- 
tisoliton solutions of the trilinear equations. Hopefully also the r-function formalism 
of the Kyoto school [16] can be extended in this direction. 
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